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Abstract 

We present long range statistical properties of a recently introduced unitary 
random matrix ensemble, whose short range correlations were found to describe a 
transition from Wigner to Poisson type as a function of a single parameter. We 
argue, by evaluating the two-level correlation function of a different solvable model, 
that the transition is perhaps a quite general feature for a class of models. Our 
analytic results compare well with numerical studies on a variety of physical systems. 
We discuss the possibility of experimentally observing signature of such transitions. 
PACS Nos. 05.45+b, 05.40+j, 71.30+h. 
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Statistical properties of eigenvalues of matrices describing a wide variety of quantum 
systems follow the universal results of random matrix models [0] . Such models were orig- 
inally conceived by Wigner to provide theoretical framework for the understanding of the 
statistics of energy levels of heavy nuclei 0, and has been developed into an important 
branch of mathematical physics through the work of Dyson and Mehta On very gen- 
eral grounds the randomness of the matrix elements, subject to any relevant symmetry 
requirements, gives rise to a model of eigenvalues repelling each other with a logarithmic 
interaction, resulting in strong correlations amongst them. The normalizability condi- 
tion of the joint probability distribution requires a confining potential for the eigenvalues, 
which can be thought of as resulting from some physical constraint (e.g. a given eigen- 
value density) 0]. For a given symmetry of the matrix, as long as the eigenvalues are 
confined well, the statistical properties of the levels in the bulk of the spectrum seem to 
be independent of the particular choice of the constraint 0, and follow universal distri- 
bution known generally as the Wigner distribution. However, it is becoming increasingly 
evident that while the level statistics of a wide variety of systems can be described very 
well by the highly correlated Wigner distribution, many of these systems show a tran- 
sition to a completely uncorrelated Poisson distribution when some relevant parameter 
is changed Such transitions in the spectral statistics might correspond to e.g. a 
chaotic-regular or a metal-insulator transition in the system. Attempts have been made 
to describe the transition in one particular statistical property, namely the nearest neigh- 
bor spacing distribution (which is sensitive to only the short range correlations between 
eigenvalues), either by purely heuristic interpolation schemes or by considering inter- 
mediate regimes where the phase space is partly chaotic and partly regular 0. These 
results differ qualitatively from the one case where the transition has been studied in 
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detail numerically, namely the case of metal to insulator transition in a disordered system 
described by the microscopic random tight binding Anderson Hamiltonian ||^. As far as 
we know, no attempt has been made to explain even heuristically the transition in other 
statistical properties, such as the number variance, Vn, or the so called A3 statistics which 
provides a quantitative measure of the long range rigidity of the spectrum, although there 
exist numerical evidences for such transitions in both chaotic and disorderd systems 



10, 11 . 



Recently we have argued that the appropriate random transfer matrix model related to 
disordered conductors belong to a new family of random matrices; a solvable model then 
predicted a very specific type of transition in the nearest neighbor spacing distribution 



as a function of a single parameter |T2[. In the present work we first show, by evaluating 
the two-level correlation function of a very different (though related) but still solvable 
model, that the nature of this transition in spectral statistics is not peculiar to the model 



considered in but is perhaps generic for a class of models with weak confining po- 
tentials. We then calculate, for the above models, two other statistical properties which 
describe the transition in the long range correlation in the bulk of the spectrum, namely 
the number variance and the A3 statistics. We show that the variance of a linear statistic 
is controlled by a single parameter in these models, and is no longer universal. Our result 
includes the theorem of Dyson and Mehta on the universality of the variance of any 
linear statistic as a special case in the appropriate limit. As an example, we give explicit 
expression for the number variance. The nature of the transition in the A3 statistics 
agrees with earlier numerical results for transfer matrices in disordered systems pUj]. 



It is well known that the Wigner distribution is satisfied in the appropriate range by 
eigenvalues of transfer matrices, Hamiltonians as well as evolution operators |]l|. Same 



is true at the other end of the transition, the Poisson distribution. Such universahty 
in the statistics results from the fact that once the symmetries of the matrix have been 
properly taken into account, most of the differences in the detail appear in the local 
density of states, which disappear when the distributions are described in units of the 
average local spacing. It is therefore quite possible that the transition from Wigner to 
Poisson distributions also has some broad general characteristics, although there may be 
more than one class of transition depending on the physical nature of the problem. It is 
therefore of interest to compare the predictions of the above mentioned solvable model 
with numerical results for different statistical properties of different systems, described 
by different kinds of matrices. We find that the transition in the A3 statistics describes 
qualitatively well the numerically obtained statistics for a variety of physical systems 
TD| , , including the energy levels of disordered systems [|rT[] . Moreover, the transition in 



the eigenvalue-spacing distribution evaluated earlier for the model also shows remarkable 
similarity to the transition in the distribution of energy levels numerically obtained in ref. 
PI mentioned above. 

We therefore conclude that in particular, the numerical results for both the nearest 
neighbor spacing distribution and the long range A3 statistics of the energy levels of 
disordered systems are well described by the solvable models, even though the models were 
originally constructed for transfer matrices. The solvable models then allow us to calculate 
and predict further consequences of such transitions in the energy level statistics. In 
particular we predict how the "correlation hole" in the Fourier transform of the absorption 
spectrum of small metallic particles should be destroyed with increasing disorder, and 
discuss the possibility of observing it experimentally. 

As mentioned earlier, the random matrix models are characterized by the confining 



potential V{x). The statistical properties of the levels can then be evaluated from the 
two level correlation function which can be obtained explicitly from a set of orthogonal 
polynomials defined with the potential as the weight factor [0]. We will call the random 



matrix ensemble introduced in ref. [0 the "g— Hermite" Unitary Ensemble, because 



the orthogonal polynomials defined by the potential are the "q-generalization" [|T^ of 

the classical Hermite polynomials that characterize the conventional Wigner or Gaussian 

Unitary Ensemble (GUE). The correlation function for the g— Hermite ensemble was found 

to depend crucially on some parameter /3 = ln(l/g) characterizing the potential: 

/3 simriu-v) 
^(^-"^ = 2;^ sinh(/j(.-.)/2) ' 

where the scaled variables u and v are such that the density K{u, u) is unity. As noted 

in ref. [1^, this reduces to the corresponding correlation function for the GUE in the 



limit /? 0. The above expression for the kernel is valid in the bulk of the spectrum 
and for < /5 < 27r^. We observe here that the nearest neighbor spacing distribution as 
a function of the parameter obtained earlier for this model is remarkably similar to the 
one obtained numerically in |]^ for the energy eigenvalues corresponding to a microscopic 
Anderson model for disordered conductors going through a metal-insulator transition. 
(The difference in the power law behavior for small spacing and the precise point where 
all different curves cross is due entirely to the fact that the numerical results are for 
orthogonal symmetry, while our model has unitary symmetry.) In order to establish that 
the above two-level kernel is not a peculiarity of the g— Hermite ensemble, we first consider 
another ensemble defined by the potential, 

oo 

\/(x;g) = ^ln[l + (l-g)g^x], (2) 

n=0 

where the eigenvalue x is from to oo as opposed to the range — oo to +oo for the 



g— Hermite ensemble. For large x and small q, this potential behaves as [Inx]^, as in the 
g— Hermite case, but for small x it behaves linearly as opposed to the quadratic dependence 
in g— Hermite. This potential was considered in as a possible model for the transfer 
matrix describing disordered conductors. The orthogonal polynomials for this potential 
are the g— Laguerre polynomials, a generalization of the classical Laguerre polynomials 



15| ; we therefore call it the "g— Laguerre" model. From the asymptotic properties of the 



g— Laguerre polynomials |T4|, we obtain the two- level kernel in the bulk of the spectrum 



as a function of /3 = ln(l/g) ^ 1 in the limit where the number of eigenvalues N — > oo: 

K{x,y) 

smf-lnxlcosf-Wj-y ^m f -In,/ 1 cos f -Inx 



x~y 



(3) 

In order to compare with the Wigner distribution for which the average spacing between 
adjacent levels is unity, it is necessary to use a transformation of variables x = e^^^, y = 
e^'^", such that the density in the new variable is uniform and unity. In this new variable, 
the kernel becomes precisely the same as given in Eq.(l) for g— Hermite ensembles. It 
therefore follows that the statistical properties of the levels in the two models are identical 
in the bulk. Note that for (3 —>■ 0, the two models reduce to classical Hermite and Laguerre 
models, which are known to have identical kernels in the bulk [Q. (We shall not consider 
here the interesting "edge effects" in the Laguerre ensemble |16|.) Thus the statistical 
properties of the levels in the bulk of the spectrum for large are insensitive to the details 
of the model as long as some general features of the model (weak [Inx]^ confinement of 
the large eigenvalues) remain the same. 

We now evaluate the number variance and the A3 statistics for the g— Hermite model, 
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expecting that these results will be valid for at least a class of similar models. The variance 
= — of the number of eigenvalues n in an interval (— s/2, s/2) can be expressed 



as 



VJs) 



s/2 rs/2 2 1 rnsA"'? 

dx dy[5{x-y)-Y{x-y)] = - dk [l-h{k)], (4) 

-s/2 J-s/2 TT Jo fc^ 



where the two- level form factor b{k) is the Fourier transform of the two- level cluster 
function Y{x) = [i^(a;)]^. For the kernel given by (1), we obtain 

47r6(fc) = (|A;| + 27r)coth[(|A;| + 27i)n/p] + i\k\- 27r)coth[(|A;| - 27r)7r//3] - 2|A;|coth(|A;|7r//3). 

(5) 

Note that b{k) reduces to the exact GUE result in the limit (3 —>■ 0. To an excellent 
approximation, b{k) contributes negligibly to the integral in Eq. (4) for k > 2ti. We find 



Vn{s) = 1 - p S + 4 [CI{2'KS) + 1] + — / dk 



where Ci is the cosine integral, and p defined as 



coth 



fcvr 



13s 
kn 



(6) 



p = coth 



3 2vr2 



(7) 



goes to 1(0) for (3 — > 0(oo). We can obtain an explicit expression for the number variance 
by approximating the factor [coth(|j) — under the integral by k7r/3(3s up to A; = S/^sp/vr 



and by p for larger k. The result 



K(s) = [i-p]s + 



P 



sin(3/5ps/7r) 
SPps/iT 



+ 



P 



In 



agrees with the numerical evaluation of the integral up to a small f3 dependent constant for 
large s as well as corrections of order 1/s. Eq.(8) clearly exhibits the crucial dependence 
of the number variance on the parameter (3. As expected, (3 = corresponds to the exact 



GUE result, with a logarithmic dependence on s; increasing /5 corresponds to a transition 



towards a Poisson result, which is hnear in s. Note that the variance of an arbitrary hnear 
statistic / = J2nf{^n), of which \4 is a special case, will depend on the form factor (5) 
and is clearly no longer universal, 
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Var(/) = -/ dk[l-h{k)]\f{k)\\ (9) 

TT JO 

where f{k) is the Fourier transform of f{x). This result generalizes the Dyson- Mehta 
theorem[0 on the variance of arbitrary linear statistics and reduces to it when /5 ^ 0. 

The A3 statistics is a measure of the size of fluctuations of a given level sequence 
against a best straight line fit for that level sequence. If s{x) is the staircase function for 
a level sequence (in the variable where the density is uniform) in a given interval, then 
one defines a variance A3 =< minA,B[(s(u) — Au — B)"^] >, where <> denotes an average 
over an ensemble of level sequences. This can be expressed in terms of the two-level kernel 
K{u, v) 0. For the g— Hermite ensemble we obtain the small and large limits explicitly: 



5 

ln(27rs) + 7 ~ ^ 



2„2 



7271 

for /3 <^ which reduces to the GUE result for /3 = 0, and 

P 2 



+ 0[{(3s)') (10) 



A3(g,/3) _ 1 
s 15 



+ C + 0(e-'^^), (11) 



exp (^) - 1 

for f3 ^ 1/s, which reduces to the Poisson result for /5 — 00. Here C is a numerical 
constant independent of s and very weakly dependent on f3. Note that since the ratio A3/S 
is finite, A3 always has a linear dependence for large enough s, with a slope that increases 
with increasing f3, approaching the Poisson limit for /3 — 00. Figure 1 shows the complete 
solution (obtained from numerical evaluation of the integrals involved) for various values 
of (3. The deviation from the GUE result as a function of some parameter is qualitatively 



similar to the deviations seen numerically in both the transfer matrix [10| and energy 



eigenvalues [|ri[] corresponding to the tight binding Anderson Hamiltonian for disordered 
conductors, as well as for the eigenvalues of the evolution operator corresponding to the 
Fermi-acceleration model [^]. 

We find that the statistical properties of the g— Hermite (or the g— Laguerre, in the 
bulk) ensembles are similar to those of a wide variety of quantum systems, including 
energy levels of disordered conductors, studied numerically. Although we are not able to 
derive these ensembles from microscopic Hamiltonians at present, it is useful to explore 
specific spectroscopic signatures of the ensembles so that their relevance to given physical 
systems can be tested experimentally. The obvious problem of looking for evidence of these 
transitions is the difficulty to extract a "stick spectrum" from experimental data where a 
large number of levels are usually lost either in the noise or in unresolved bands. One way 
to avoid this problem is to take direct Fourier transform (FT) of the raw experimental 
data; the ensemble average of the square of the FT will show a "correlation hole" if the 



spectrum is chaotic, i.e. if it has Wigner distribution [0. The size of the correlation hole 
is proportional to [1 — b{k)], where b{k) is the two level form factor mentioned before. We 
plot this function for the g— Hermite ensemble obtained from eq. (5) for various values 
of (3 in figure 2. The depth of the hole decreases from the Wigner result in a specific 
way towards the Poisson limit of no correlation hole {b{k) = 0). Thus e.g. microwave 
experiments on ensembles of small metallic particles (of roughly equal size) |jl8| at various 
disorders might reveal this behavior (increasing disorder will correspond to increasing jS 
r2|). Of course one needs sufficiently low temperature and small system size so that the 
energy levels are not broadened into a continuum. For a metallic system, a crude estimate 
within a simple electron gas model suggests that for 50 nanometer size particles one might 
expect to observe the effect at a temperature in the milikelvin range. 
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Figure captions: 

Figure 1: A3 statistics as a function of the length s of a level sequence for different 
values of (3 for the g— Hermite model. The (3 = curve coincides with the GUE result 
while the (3 = 00 line coincides with the Poisson result. 

Figure 2: The function [1 — &(A;)] showing the change in the correlation hole with 
increasing (3. Interpreting k as time and with the dashed line replacing an ideal vertical 
drop at very small time for finite number of levels n, this should correspond to the square 
of the Fourier transform of an experimentally obtained spectrum for a system of small 



disorderd particles. Increasing (3 will correspond to increasing disorder . 
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